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It is argued that two-particle collisions of relativistic particles "at a distance", irrespective of their position 
I in the configuration space, generate uniform distribution of particles in Lorentz invariant phase space. 

. In a recent paper by Meres at al. [1] a new numerical procedure for generation of relativistic 
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particles with the given total energy-momentum is proposed. The procedure makes use of two- 
particle collisions between particles considered in momentum space only, with no definite positions 
in configuration space, and seems to produce uniform distribution of particles in Lorentz invariant 



pi • phase space (LIPS) when applied repeatedly to a system with non-uniformly distributed particles at 

D ; 

r~| \ the beginning. The paper includes a sketch of an argument explaining why it is so. In this note a 



detailed exposition of the argument is given. 

In the procedure developed in [1], collisions are linked into sequences with two collisions per particle 
and the pairing of particles in each sequence is governed by a single parameter whose value is chosen 
at random. While this can be advantageous in numerical calculations, it is technically simpler and 
physically more instructive to consider a gas of particles colliding completely at random. Our aim 
is to show that the equilibrium distribution of particles in such a gas is uniform in LIPS, so that 
the distribution necessarily becomes uniform asymptotically^ for the number of collisions approaching 
infinity. If this is the case, the procedure with linked collisions, which uses essentially the same 
mechanism, will yield an asymptotically uniform distribution, too. 

Consider a gas of nonidentical relativistic particles in which there are perpetually going on 
elastic 2-particle collisions with a given angular distribution of the outgoing particles in cms. The 
particles are free between collisions and collide "at a distance", irrespective of their actual position in 
the configuration space. To describe the evolution of the gas we divide it into steps with one collision 
per step, and pick the pair of particles entering the collision in any given step at random. Then we 
choose the directions in which the particles are coming out of the collision in accord with their angular 
distribution in cms. Both choices are obviously independent of the choices made in previous steps, as 
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well as of the sequence number of the actual step. Thus, the evolution of the gas is a homogeneous 
Markov process. 

There is a vast body of literature on Markov processes, including monographs devoted entirely 
to them (see, for example, [2]) and books on probability theory where their properties arc discussed 
in detail. For a newcomer, the best choice is perhaps to read an exposition on the subject in some 
lecture notes available online, as are those cited in [3]. 

A homogeneous Markov process is described by the transfer matrix Pij, whose element {ij) is the 
probability that the system passes in one step from the state j to the state i. Denote the probability 
that the system will be found in the nth step in the state i by Pn{i)- The probabilities in two 
subsequent steps are related by the formula 

Pn+iii) = ^PijPn{j)- 

3 

In statistical mechanics, one calls the states i, j, ... "microstates", and the states defined by the prob- 
ability distribution "ensemble states'' or "macrostates". Equilibrium is a macrostatc into which 
the system gets after infinitely many steps from the initial macrostate with an arbitrary probability 
distribution; in particular, an arbitrary microstate (which can be regarded as a macrostate with one 
probability equal to 1 and all others equal to 0). Thus, for the probability distribution in equilibrium 
we have 

p{i) = y,ViMj)^ ^= 1™ (1) 

^ — ' n— >-oo 
j 

On the other hand, equilibrium can be defined also as a macrostate whose probability distribution 
stays constant during the evolution of the system. Thus, the probability distribution in equilibrium 
should satisfy 

p{i) = Y,P^oP{j). (2) 
j 

In the first lecture notes cited in ref. [3], the interplay between equations (1) and (2) is demonstrated 
on a simple, but no less illuminating, example of two-level system. 

According to equation (2), the vector p is an eigenvector of the matrix P with the eigenvalue 1. 
One can easily see that there necessarily exists some vector with this property, the reason being that 
the matrix P obeys the normalization condition 

Y,Pi3 = ^- (3) 

i 
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Indeed, if we sum both sides of equation (2) over i and use equation (3), we arrive at identity; thus, 
the determinant of the system is zero and the system has at least one nontrivial solution (modulo 
rescaling). However, we would also like to know that the solution is unique and nonnegative, and that 
it can be expressed in the form (1). It turns out that all this is indeed the case provided the matrix 
P is irreducible; that is, either the matrix itself or some finite power of it is positive (has all entries 
positive). In less formal terms, the solution has desirable properties if we are able to pass from any 
state to any other state in a finite number of steps, so that we do not eventually become trapped in 
some subspace of the state space. The proof is based on a well-known theorem from the theory of 
nonnegative matrices, called Perron-Probenius theorem. More details can be found in [4]. 

The microstates of the gas are numbered by the A^'-tuples of particle momenta (pi, p2, . . . , Pn)- 
Thus, there are continuously infinitely many of them. However, we can pass to the discretized theory, 
dividing the momentum space into elementary cells of finite size. This is actually done in numerical 
calculations. In the discretized theory, the number of microstates is countable, and if we take into ac- 
count that the probability distribution we are interested in includes a 5-function on total 4-momentum 
(which is replaced by l/[4-dimensional volume of the order (size of the cell)^] x Kronecker delta after 
discretization), the number of microstates becomes finite. As a result, previous considerations can be 
applied without modifications to a gas of colliding particles as a special case. 

The basic building block of the theory is the expression for the distribution of outgoing particles in 
a two-particle collision. Let us compute this distribution in the laboratory frame in case it is isotropic 
in cms. The quantities referring to cms will be denoted by the index 0. The probability that the 
particle 1 comes out of the collision in cms in the solid angle d^io is 

,„ d^io 
dP = — — . 
Ait 

Denote the magnitude of the momenta of ingoing particles in cms by ttq, the momenta of outgoing 
particles in cms by (pio,P2o)) the total energy of ingoing particles in cms by cq, and the energies 
of outgoing particles in cms by (£'10,-^20 )• The energies are defined in terms of momenta as £'10 = 



\/Pio + ^\ 9-^*^ -^20 = \/P20 + "^2' where m\ and m2 are the masses of the colliding particles. The 
expression for dP can be rewritten as 



dP = - — ^5{pw - 7ro)d^pio = -; — ^5{pu) - 7ro)5(pio + V>2Q)d^Pwd^P20, 



and if we use 



S{Eio + E20 - eo) = S{pio - TTo), 
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we find 

dP = V-glO + 1/^20 ^^^^^ ^ _ g^)^(p^Q + P20)dV0d^P20. 

47rpio 

The product of (5-functions can be written as (5^(pio + P20 ~ Qo)^ where p^Q = (-EiOiPio) ^^id P20 ~ 
(^20)P2o) are the 4-momenta of outgoing particles in cms and Qq = (eo,0) is the total 4-momentum 
of ingoing particles in cms. Prom the conservation laws we also have l/£'io + l/-£^20 = eo/(-£'io-S2o) 
and pio = TTo, hence 

dP = (Pio +P20 - go) p p ■ (4) 

47r7ro -C'io-C'20 

The 4-dimensional (5-function as well as the element of LIPS by which it is multiplied are relativistic 

invariants, so that they can be both carried over into the laboratory frame simply by wiping off the 

index 0. The factor in front of the 5-function must be rewritten as an invariant, too, which means that 

we must express ttq in terms of eo and subsequently replace eo by (invariant mass). By squaring 

the equation eo = Eiq + E20, separating the term 2E10E20 and squaring the resulting equation again 

we obtain 

^0 - 2eo(£'io + -E20) + (-^10 ~ -^20)^ = 0' 
and if we insert here the expressions for £'10 and £^20 with pio = p20 = ttq and solve for ttq, we find 

TTO _ 1 A _ 2M2 «4y/2 

where = ml + 771-2 and fj? = m\ — m^. Inserting this in equation (4) and passing from cms to 
laboratory frame we finally obtain 

1 / 2M2 un-1/2 



The factor in front of the (5-function in the expression for dP serves as a normalization factor, making 
the total probability equal to 1. The correct normalization of dP is guaranteed by the way the 
expression has been derived, and can be easily verified in the nonrelativistic limit. 

If the distribution of outgoing particles in cms is anisotropic, we must insert a factor /(cos ^10) 
into the expression for dP, where ^10 is the angle between the initial and final directions of motion of 
the particle 1 in cms and / is a nonnegative function satisfying 

1 

f{x)dx = 2. 



/■ 
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The argument of / can be expressed in terms of and q\ . pi, where qi and pi are the initial and 
final 4-momenta of the particle 1. In particular, for mi = m2 = m we have 

cos 6^10 = — 5 — 

— 4m^ 

Consider now a gas of N particles with momenta (qi, q2, ■ ■ ■ , Qat), and suppose that at the given 
moment there occurs a two-particle collision, with all pairs of particles colliding with equal proba- 
bility. We are interested in the probability of finding the gas after the collision in an infinitesimal 
neighborhood of the point (pi,P2, ■ ■ ■ , Pat) of the iV-particle momentum space. This can be written 
as a sum of contributions of all pairs of particles, 

dP = JfYl '^^'^ n <^(P'= - ^c)d^Pc, (6) 

a<b c^a,b 

where M = N{N — l)/2 is the number of pairs and dPab is the probability that after the particles 
a and b have collided, their momenta will be found in an infinitesimal neighborhood of the point 
(pa, Pb) of the two- particle momentum space. Suppose for simplicity that the distribution of outgoing 
particles in 2-particle collisions is isotropic in cms. Then the probability dPab is given by equation (5) 
rewritten for particles a and b, 

V Qtb %.b J 

where Pab = Pa+ Pb, Qab = Qa + Qb, M^b = "^a + "^6 and /x^^ = ml - ml. Taking into account the 
(5- functions in equation (6), we can rewrite the 5- function in equation (7) as S^{p — q), where q'^ and 
arc the total 4-momcnta of in- and outgoing particles. We can also transform the noninvariant 
3-dimensional (5-functions in (6) into invariant ones by multiplying them by the energies of the ingoing 
particles ec- In this way we arrive at the formula 

dP = F5\p~q)ll^, (8) 

with the distribution function F defined as 

The probability distribution dP, regarded as a function of initial and final momenta of the particles, 
is the transfer matrix for our problem. It is continuous, but can be discretized by the procedure 
mentioned earlier. 
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The gas in equilibrium is described by the probabihty distribution 

dp = f6\p-P)ll^, (10) 

where p^^ is the total 4-momentum of the gas, is the value assigned to p^ and / is the distribution 
function in equilibrium depending on the momenta (pi, P2, ■ ■ ■ ,Pn) of the particles forming the gas. 
Denote the A^-tuple of momenta (pi,P2, • • • ,PAr) by p. For the probability distribution dp we have 
the continuous version of equation (2), 

dp{p) = J dP{p,q)dp{q), 

where we integrate over the A/^-tuple of momenta q. By inserting here from equations (6) and (10) we 
obtain 

fip)6\p -P)= [ Fip, q)5\p - q)f{q)5\q - P) [] 

and if we replace q^ in the argument of the second (5-function on the right hand side by p^, using the 
fact that the value of q^^ is fixed by the first (5-function, we find 

f{p)= [ F{p,q)f{q)6^ip-q)ll^. (H) 

After having obtained an integral equation for the distribution function /, our next goal is to 
prove that it is solved by constant /. Thus, we want to show that F satisfies the equation 

fF{p,q)6\p-q)l[^ = l. (12) 

Note that a similar equation with the integration over the A/^-tuple of momenta p, 

lF{p,q)6\p-q)l[^ = l, 

is just the normalization condition for F. In fact, it is a continuous version of the normalization "by 
columns" introduced in equation (3). Our claim is that the transfer matrix we have constructed is 
normalized also "by rows". (Strictly speaking, the matrices in the two normalization conditions are not 

the same, since they have different energies in the denominator. [The elementary volumes referring 
to different momenta do not matter, since they can be always chosen to have a fixed size throughout 
the momentum space.] However, we can view both matrices as different representations of a single 
matrix.) 

To prove equation (12), we first transform it into an equation for two-particle probability distri- 
bution. For that purpose, we pick an arbitrary term (a, b) in the sum over the pairs of particles in 
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(9), rewrite the energy-momentum conservation law so that only the 4-momenta of particles a and 
b remain in it, and eliminate the 3-dimensional (5- functions by integrating over all q's other than 
(qa, q;,). In this way we obtain 

/ 7^ 1 — + — ^ iP^2 - Qu) = 1- (13) 

This can be viewed as normalization of the two-particle transfer matrix "by columns". To show that the 
equation is satisfied, it suffices to replace in the expression in front of the (J-function by p^2) rename 
(pi,P2) to (qi,q2) and vice versa, and notice that the resulting equation is just the normalization 
condition for the probability dP defined in equation (5). Thus, the normalization of the two-particle 
transfer matrix "by columns" follows immediately from its normalization "by rows". 

Previous reasoning generalizes trivially to anisotropic collisions. In this case, the integral in equa- 
tion (13) contains an additional factor /(cos^io), but the equation is still satisfied because the in- and 
outgoing momenta enter the expression for symmetrically. This is seen in the formula for cos ^lo 
cited before, and can be regarded as a consequence of the trivial fact that the angles by which the 
vector pio is defiected from the vector qio and vice versa are identical. Thus, equation (12) is valid 
no matter what the distribution of outgoing particles in cms. 

We have shown that the condition of equilibrium is solved by the probability distribution that is 
uniform in LIPS, 

dp = CpS\p-P)l[^, (14) 

where Cp is a normalization constant depending only on the 4- momentum of the gas P^. Thus, the 
particles forming the gas, if in equilibrium, can be distributed uniformly in LIPS. However, we would 
like to be sure that they are distributed uniformly in LIPS. In the discretized theory this is guaranteed 
by the Perron-Probenius theorem, provided the gas can pass from any initial state to any final state 
after a finite number of two-particle collisions. In fact, the final state does not need to be varied, since 
the desired transition for any initial and final state can be combined from one direct and one inverse 
transition to some reference state. Wc will not attempt to prove that the interpolating sequence of 
collisions exists in general, but we will provide an illustrative example showing that such sequence 
exists for a pair of states that are in a sense maximally distant from each other, and that its length 
is of order N. 

Consider a many-particle gas with the total energy E and zero total momentum, and suppose for 
simplicity that the particles are massless. We are interested in the transitions into the reference state 
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that is maximally randomized, with the energy of all particles equal to = E/N and the momenta 
filling uniformly a sphere in the momentum space. If the gas is to evolve into this state, the least 
favorable case seems to be if it is in the state with a certain particle A having energy 8q = E/2 and 
momentum oriented in one direction, and all remaining particles having energy e = £q/{N — 1) = 6^/2 
and momenta oriented in the opposite direction (fig. 1 left). Consider a sequence of A^ — 1 collisions in 




Fig. 1: Randomization of a state with half the total energy carried by a single particle 

which the particle A pairs subsequently with all remaining particles, knocking them alternately to one 
and the other side in a fixed plane and setting their energy to e^. The energy-momentum conservation 
laws written for the i^th collision are (fig. 1 right) 

£' + er = £ + e, 8' cosa + ercos f3 = £ — ecosO, sin a + e,- sin/3 = e sin 0, (15) 

where £ and £' is the initial and final energy of the particle A, a is the scattering angle of the particle 
A, /3 is the angle between the initial momentum of the particle A and the final momentum of the 
particle B that is colliding with A, and 6 is the cumulative scattering angle of the particle A (the 
angle between the initial momentum of A in the z^th collision and the momentum A had before the 
first collision). In all scatterings but a couple of last ones the energy £ is much greater than both 
e and e' (the final energy of the particle B) = 2e. This causes a to be small and /3 to be close in 
absolute value to tt/2. If we choose the first particle B to be declined from its original direction to the 
right of the particle A, (3 will be of the form (3 = =F7r + 5, where 8 is small and the upper and lower 
sign correspond to odd and even collisions respectively. (In the figure, an odd collision is depicted.) 
As for the third angle appearing in the problem, for a sequence of collisions under consideration, with 
an altering direction of motion of particles B after collision, is always less in absolute value than 
the actual a, hence it is small, too. 

We want to solve equations (15) in the leading order in the small parameter e/£. This is a simple 
exercise, however, at one point we must be cautious: we can put = 2e almost everywhere, but in 
the term fees a in the second equation (15), in which appears through £' ^ we must use the exact 
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formula = 2(1 — \/N)e. Prom the first equation (15) we obtain £' = £ — e, which yields 

£' = So — ve = fie, (16) 

where fi is the number of collisions which remain till the end. For our purposes we can put ji = N — u. 
(To obtain correct values of £' at the last stage of the process, we must pass loiJb = N— 1 — v and 
account for the cumulative effect of the small correction to e^.) From the remaining two equations we 
find 

« = ±-, 5 = ±{--^)- (17) 

To compute 9, wc must sum q's for all collisions up to that with the sequence number u — 1. For odd 
collisions it holds 

1 



/, 2 2 2 
^odd = — 7 - — + • • • + 



Mm + 1) 

and for even collisions we have 



A^-1 N-2 II fi + 1 

N-l 

1 1 

+ ... + 



(A^-l)(iV-2) 



+ 



dx__l 1 ^ _l 1_ 



112 1111 
+ TT+ r = + — = - + —. 



fi + l N 11 + 1 11 + 1 N 11 N 



The momentum acquired by the particle B after it collided with the particle A is deflected from the 
direction perpendicular to the original momentum of the particle A by the angle A = 5 + 9. After 
inserting here for 5 and 9 we obtain 

2 

Aodd = 0, Aeuen = (18) 

The approximate theory developed here holds for a majority of the particles of the gas. After 
the collisions with the particle A are completed, these particles split into two equally large groups, 
one containing particles knocked to one side of the particle A and the other containing particles 
knocked to the other side. From equation (18) we can see that the particles of the first group move 
perpendicularly to the original direction of motion of the particle A, while the particles of the second 
group are slightly deflected to that side to which the particle A was heading at the beginning. A 
small third group consists of particles that collided last and move in various directions in such a way 
that their total momentum balances the momentum arising from the deflection of the particles of the 
second group. 



Collisions with the particle A represent just the first step of a three-step procedure that brings 
the gas to the reference state. The next step is to modify the direction of motion of the particles 
so that they combine into pairs with exactly opposite momenta; or in the discretized theory, with 
momenta that lie in exactly opposite cells of the momentum space. If the size of the cells Ap is 
greater than ApQ = Aeven^r = 4e/A/', the task is restricted to the particles of the third group and 
can be presumably completed in a number of steps of order of the number of particles belonging to 
that group, which is negligible in comparison to N. In the opposite case we must take care of the 
first two groups of particles, too. We must slightly rotate, say, the direction of motion of the particles 
of the second group so that it will become opposite to the direction of motion of the particles of the 
first group. This is easily accomplished with one collision per particle; we just have to think of it in 
advance and pick a slightly different final state in the collisions of this group of particles with the 
particle A. A straightforward computation yields that the final energy e' of the particles of the second 
group must differ from the energy by the quantity 

Ae' = ^^""^'^ (19) 

A/"(l -|-sin^ + cos<;ii)' ^ ^ 

where (f) is the angle by which the momentum of the particle of the third group, paired with the given 
particle of the second group, is deflected from the original momentum of the particle A. Note that the 
expression for Ae' diverges for = vr, therefore the procedure does not work with particles that have 
not collided with the particle A yet. To find how many collisions arc necessary, notice that we must 
straighten the momenta of half the number of particles, and need one collision for one straightening 
(perhaps a little more for the particles of the third group, but their number is negligible anyway). 
Thus, the number of colhsions needed at this stage is N/2. 

The third and final step of the randomization procedure is to distribute the particles uniformly on 
a sphere in the momentum space with the radius e^. This is done simply by colliding particles with 
opposite momenta. The number of collisions is again N/2, which makes the total number, in the less 
favorable case when the cells in the momentum space have a size less than Apo, equal to 2N. 

An isolated gas whose particles collide with each other wanders the energy hypersurface of the 
A'^-particle momentum space. If the hypersurface consists of domains that are interconnected but not 
connected with each other, in the sense that points of one domain cannot be reached from another 
domain, there can exist a macrostate that stays unchanged once it was established, but the gas will 
not evolve towards it from other macrostates. The question is whether the uniform distribution in 
LIPS defined in equation (14) does not represent such a "virtual" equilibrium state. The example with 
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randomization of the gas with one highly energetic particle disputes this possibility since it shows that 
even faraway points of the momentum space can be accessible to each other. 

For a given homogeneous Markov process, one can calculate the equilibrium probability distribu- 
tion numerically by starting with an arbitrary state Iq, booking the state after every S> 1 steps 
and computing the relative frequency of every state after v ^ steps. The probability distribution 
Pappii) obtained in this way is arguably a good approximation to p{i). To see that, consider the first 
sequence of n steps. The probability distribution p^{i) at the end of the sequence is given by a cut-off 
version of equation (1) with po{i) = Su^, 

3 

For large enough ^ the information about the initial state is effectively wiped out during the sequence, 
so that the effect of the cut-off is small and p^{i) is close to p{i). Thus, we could find Pappii) simply by 
calculating Pn{i) numerically; that is, by repeating the sequence with the same initial state sufficiently 
many times and computing the relative frequencies of the final states. However, if we picked by chance 
an initial state for which the system would wander through marginal regions of the state space for 
a long time (this can presumably happen for initial states with small probability), the resulting 
probability distribution would be distorted. To be on the safe side, we can use the procedure outlined 
at the beginning, with the final state of one sequence serving as the initial state of another sequence. 
This can be described by the formula 

Pappii) = ^V\fpapp{3). 

i 

The same formula is valid, as a trivial consequence of equation (2), for p{i). Nevertheless, Papp{i) is 
still an approximation, since the number of events used in its computation is finite. For a relativistic 
gas, numerical calculations following this recipe support our claim that the equilibrium distribution 
of particles is uniform in LIPS [5]. 

To complete the discussion, consider a gas with linked collisions introduced in [1]. Suppose the two- 
particle collisions come in sequences consisting of N collisions, with the sequence number of colliding 
particles, regarded as a periodic variable with the period N , raised by one in each consecutive collision. 
Furthermore, suppose that the difference between the sequence numbers d, which stays constant 
throughout each sequence, assumes all values 1, 2, . . ., A'' — 1 with the same probability. The collisions 
fall into three classes (fig. 2): colhsions No. 1 to D = min(cZ, N — d) ui which both incoming particles 
are in the initial state (class a), collisions No. D + \ io N — D (absent for N even and D = N/2) 
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in which one incoming particle is in the initial state 




and the other incoming particle is in the intermediate 



state (class /3), and collisions No. N — D -\- 1 to in 



which both incoming particle are in the intermediate 
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state (class 7). The two-particle probability distribu- 
tions in the respective classes are dP^") = dP{qq,ll), 
dpW = dP{ql,lp) and dP^ = dP{ll,pp), where q, I 



and p are initial, intermediate and final momenta re- 



spectively. The A^-particle probability distribution for 



5 



one sequence of collisions is 



dP = 



1 




(20) 



A^- 1 



d, 



Fig. 2: Pairing of particles in collisions 



where the contribution of collisions with the given shift 



in the sequence number can be written as a product of two-particle probability distributions inte- 
grated over the intermediate momenta, 



The A^-particle distribution is again of the form (8), and is again normalized "by columns" thanks to 
the fact that the two-particle distribution is normalized "by columns". (Equation (12) for the distri- 
bution function reduces again to the identity (13).) Thus, the linking of collisions does not affect our 
conclusion that the uniform distribution in LIPS obeys the condition of equilibrium. The integra- 
tions in (21) remove the better part of the (^-functions in the two-particle probability distributions, 
leaving us with a much greater freedom in choosing the final momenta than we had in the case with 
independent collisions. Because of that it could be expected that in a many-particle gas the particles 
will need only 0(1) sequences of collisions to get from any state to any other state. This, if true, 
would make the total number of intermediate collisions of order N, just as we have established in the 
example with randomization. However, to prove that seems to be no less difficult for linked collisions 
than for independent ones. On the other hand, if the particles could indeed pass from any state to 
any other state within 0{N) independent collisions, they could of course do that also within 0{N) 
sequences of collisions. As a result, if the true (not "virtual") equilibrium distribution is uniform in a 
gas with independent collisions, it is necessarily uniform also in a gas with linked collisions. 
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